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Abstract. We discuss the implementation and results of a recently developed microscopic method for calcu- 
lating ion-ion interaction potentials and fusion cross-sections. The method uses the TDHF evolution to obtain the 
instantaneous many-body collective state using a density constraint. The ion-ion potential as well as the coor- 
dinate dependent mass are calculated from these states. The method fully accounts for the dynamical processes 
present in the TDHF time-evolution and provides a parameter-free way of calculating fusion cross-sections. 



1 Introduction 

The study of internuclear potentials for heavy-ion colli- 
sions is of fundamental importance for calculating fusion 
cross-sections, and for studying the formation of super- 
heavy elements and nuclei far from stability. While asymp- 
totically such potentials are determined from Coulomb 
and centrifugal interactions, the short distance behavior 
strongly depends on the nuclear surface properties and the 
readjustments of the combined nuclear system, resulting 
in potential pockets, which determine the characteristics 
of the compound nuclear system. 

Among the various approaches for calculating ion-ion 
potentials are: 1) Phenomenological models such as the 
Bass model 11], the proximity potential |2|, and poten- 
tials obtained via the double-folding method |3,4|. Some 
of these potentials have been fitted to experimental fusion 
barrier heights and have been remarkably successful in 
describing scattering data. 2) Semi-microscopic and full 
microscopic calculations such as the macroscopic - mi- 
croscopic method ||5]|6l, the asymmetric two-center shell- 
model |7|, constrained Hartree-Fock (CHF) with a con- 
straint on the quadrupole moment or some other defini- 
tion of the internuclear distance |l8]|9l, and other mean-field 
based calculations 1 101I11II12I . 

One common physical assumption used in many of the 
semi-microscopic calculations is the use of the frozen den- 
sity or the sudden approximation. As the name suggests, 
in this approximation the nuclear densities are unchanged 
during the computation of the ion-ion potential as a func- 
tion of the internuclear distance. On the other hand, the 
microscopic calculations follow a minimum energy path 
and allow for the rearrangement of the nuclear densities 
as the relevant collective parameter changes. In this paper, 
we shall call this the static adiabatic approximation since 
a real adiabatic calculation would involve a fully dynami- 
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cal calculation, thus also including the effects of dynamical 
reaiTangements. 

Recently, we have developed a new method to ex- 
tract ion-ion interaction potentials directly from the TDHF 
time-evolution of the nuclear system. In the DC-TDHF ap- 
proach 1 13 1 the TDHF time-evolution takes place with no 
restrictions. At certain times during the evolution the in- 
stantaneous density is used to perform a static Hartree- 
Fock minimization while holding the neutron and pro- 
ton densities constrained to be the corresponding instanta- 
neous TDHF densities. In essence, this provides us with the 
TDHF dynamical path in relation to the multi-dimensional 
static energy surface of the combined nuclear system. In 
this approach there is no need to introduce constraining op- 
erators which assume that the collective motion is confined 
to the constrained phase space. In short, we have a self- 
organizing system which selects its evolutionary path by 
itself following the microscopic dynamics. Some of the ef- 
fects naturally included in the DC-TDHF calculations are: 
neck formation, mass exchange, internal excitations, defor- 
mation effects to all order, as well as the effect of nuclear 
alignment for deformed systems. 



2 DC-TDHF Method 

2.1 Formalism 

The concept of using density as a constraint for calculat- 
ing collective states from TDHF time-evolution was first 
introduced in Ref. [14J, and used in calculating collective 
energy surfaces in connection with nuclear molecular res- 
onances in Ref. flSl. 

In this approach we assume that a collective state is 
characterized only by density p, and cuiTent j. This state 
can be constructed by solving the static Hartree-Fock equa- 
tions 

< 0pj|ajap//|0pj >= , (1) 
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subject to constraints on density and current 

< ^pjIpWI^pj > =p(r,f) 
<0pj|j(r)|0pj>=j(r,f). 

Choosing p(r, f) and j(r, t) to be the instantaneous TDHF 
density and current results in the lowest energy collec- 
tive state corresponding to the instantaneous TDHF state 
1 0(f) >, with the corresponding energy 



(2) 



This collective energy differs from the conserved TDHF 
energy only by the amount of internal excitation present in 
the TDHF state, namely 



E*{t) - EjDHF - Ecoll(t) ■ 



(3) 



However, in practical calculations the constraint on the cur- 
rent is difficult to implement but we can define instead a 
static adiabatic collective state |0p > subject to the con- 
straints 

< 0p\pir)\0p > =p(r,f) 

< '^plyWI^p > = . 

In terms of this state one can write the collective energy as 



-■coU 



Eki„(p(t),j(t)) + Eocipir, 0) 



(4) 



where the density-constrained energy E^c^ and the collec- 
tive kinetic energy Eki„ are defined as 



Edc = < (Pp\H\0p > 
2m 



j(tf/p(t) 



From Eq. |4]is is clear that the density-constrained energy 
Edc plays the role of a collective potential. In fact this is 
exactly the case except for the fact that it contains the bind- 
ing energies of the two colliding nuclei. One can thus de- 
fine the ion-ion potential as 



V^EDc(p(r,t))-EA 



(5) 



where Ea, and Ea^ are the binding energies of two nu- 
clei obtained from a static Hartree-Fock calculation with 
the same effective interaction. For describing a collision of 
two nuclei one can label the above potential with ion-ion 
separation distance R{t) obtained during the TDHF time- 
evolution. This ion-ion potential V{R) is asymptotically 
correct since at large initial separations it exactly repro- 
duces VcouiombiRma.x)- In addition to the ion-ion potential 
it is also possible to obtain coordinate dependent mass pa- 
rameters. One can compute the "effective mass" M{R) us- 
ing the conservation of energy 



MiR) = 



2[£e.m. - ViR)} 

R2 



(6) 



where the collective velocity R is directly obtained from 
the TDHF evolution and the potential V{R) from the den- 
sity constraint calculations. 



2.2 Calculation of R 

In practice, TDHF runs are initialized with energies above 
the Coulomb barrier at some large but finite separation. 
The two ions are boosted with velocities obtained by as- 
suming that the two nuclei arrive at this initial separation 
on a Coulomb trajectory. Initially the nuclei are placed 
such that the center of mass is located at x = j/ = z = 0, 
and the x - z plane represents the collision plane. During 
the TDHF dynamics the ion-ion separation distance is ob- 
tained by constructing a dividing plane between the two 
centers and calculating the center of the densities on the 
left and right halves of this dividing plane. The coordinate 
R is the difference between the two centers. The dividing 
plane is determined by finding the point at which the tails 
of the two densities intersect each other along the jc-axis. 
However, this procedure starts to fail after a substantial 
overlap is reached. Instead, one can define the ion-ion sep- 
aration asR - Rq VI 620 L where 220 is the mass quadrupole 
moment for the entire system, calculated by diagonalizing 
the quadrupole tensor to obtain the quadrupole moment 
along the principal axis, and Rq is a scale factor determined 
to give the correct initial separation distance at the start of 
the calculations. Calculating R this way yields numerically 
identical results to the previous procedure until that proce- 
dure begins to fail and continues smoothly after that point. 



2.3 Fusion cross-section 

We now outline the calculation of the total fusion cross 
section using a coordinate-dependent mass M{R) and po- 
tential V(R). Starting from the quantized Hamiltonian 

H(R, P)^^ [M(Ry-2PM(Ry^-P] + V{R) . (7) 

with the momentum operator P = -ihd/dR, the total fu- 
sion cross cross-section 



L=0 



(8) 



can be obtained by calculating the potential barrier pene- 
trabilities Ti from the Schrodinger equation for the relative 
motion coordinate R using the Hamiltonian ^ with an ad- 
ditional centrifugal potential 



H(R, P) + 



h^L{L+\) 
2M(R)R^ 



MR) = . (9) 



Alternatively, instead of solving the Schrodinger equa- 
tion with coordinate dependent mass parameter M{R) for 
the heavy-ion potential V(R), we can instead use the 
constant reduced mass ju and transfer the coordinate- 
dependence of the mass to a scaled potential U(R) using 
the well known coordinate scale transformation. 



dR 



I M(R) 



dR 



(10) 
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Integration of Eq. ( fTOb yields 

R^m ^ R^r'(R). (11) 

As a result of this point transformation, both the classi- 
cal Hamilton function and the corresponding quantum me- 
chanical Hamiltonian, Eq. now assume the form 

- - p2 
H{R,P)^ — + UiR), (12) 

and the scaled heavy-ion potential is given by the expres- 
sion 

U(R) = V(R) = V{f-\R)) . (13) 

The fusion barrier penetrabil- 

ities r/,(£'c.m.) are obtained by numerical integration of the 
two-body Schrodinger equation using the incoming wave 
boundary condition (IWBC) method. IWBC assumes that 
once the minimum of the potential is reached fusion will 
occur. In practice, the Schrodinger equation is integrated 
from the potential minimum, R^m, where only an incoming 
wave is assumed, to a large asymptotic distance, where it 
is matched to incoming and outgoing Coulomb wavefunc- 
tions. The barrier penetration factor, TiiE^.m.) is the ratio 
of the incoming flux at R^in to the incoming Coulomb flux 
at large distance. Here, we implement the IWBC method 
exactly as it is formulated for the coupled-channel code 
CCFULL described in Ref. [ 16 1. This gives us a consistent 
way for calculating cross-sections at above and below the 
barrier energies. 

2.4 Fusion with alignment 

In the case of one or both of the reaction partners being de- 
formed one has to incorporate the nuclear alignment into 
the evolution of the heavy-ion collision dynamics. This is 
done in two steps |17|: a) A dynamical Coulomb align- 
ment calculation to determine the probability that a given 
nuclear orientation occurs at the distance R(to), where the 
TDHF run is initialized. The alignment generally results 
from multiple E2/E4 Coulomb excitation of the ground 
state rotational band. The distance /?(fo) is chosen such that 
the nuclei only interact via the Coulomb interaction, b) A 
TDHF calculation, starting at this finite internuclear dis- 
tance /?(?()), for a fixed initial orientation of the deformed 
nucleus. Since the experiments are usually done with un- 
polarized beams, in a full quantum mechanical calcula- 
tion one would have to average over discrete quantum me- 
chanical rotational bands. In the classical limit, this cor- 
responds to averaging over orientation angles. In the case 
of one spherical nucleus and one deformed reaction part- 
ner, the total fusion cross section is given by an integral 
over all orientation (Euler) angles, with soUd angle element 
dQ - Insinpdp 

(T{E,,^)^ ^ dQ^(r{E,,^,,Q) , (14) 

where dPjdQ represents the alignment probability and 
o"(£^c.ra., ■Q) is the fusion cross section associated with a 



particular alignment calculated using the ion-ion potential, 
V(R,/3), obtained from the TDHF collision for which the 
deformed partner is initialized with angle /3 with respect to 
the collision axis. Details for the most general case is given 
inRefs. llTlfTSl . 



3 Application 

In this section we will give selected examples of the DC- 
TDHF method for calculating fusion cross-sections. Cal- 
culations were done in a 3-D Cartesian box large enough 
to avoid any initialization or box boundary effects. We have 
used the full Skyrme force (SLy4) 1 19|. The numerical ac- 
curacy of the static binding energies and the deviation from 
the point Coulomb energy in the initial state of the collision 
dynamics is on the order of 50 - 150 keV. We have per- 
formed density constraint calculations at every 10-20 fm/c 
interval. The accuracy of the density constraint calcula- 
tions are commensurate with the accuracy of the static cal- 
culations. 
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Fig. 1. (Colour on-line) Potential barriers obtained from den- 
sity constraint TDHF calculations at three different energies. The 
three dashed curves correspond to the transformed potential using 
coordinate dependent masses. 



3.1 Spherical system 

The DC-TDHF method is expected to do best for nuclei 
that are well described by the Skyrme HE calculations. 
One such reaction is the fusion of '^0+^"**Pb system. In 
Fig. [1] we show an example of microscopic potentials for 
the '^0-(-^"^Pb system at three difference center-of-mass 
energies |20|. The dashed curves are the corresponding po- 
tentials transformed via the microscopically calculated ef- 
fective mass, M(R). We observe that all of the scaled barri- 
ers give a very good description of the fusion cross-section 
at higher energies suggesting these cross-sections are pri- 
marily determined by the barrier properties in the vicinity 
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Fig. 2. (Colour on-line) Total fusion cross section as a function 
of cm. energy using the potentials of Fig. [T] 



of the barrier peak, whereas for the extreme sub-barrier 
cross-sections are influenced by what happens in the in- 
ner part of the barrier and here the dynamics and conse- 
quently the coordinate dependent mass becomes very im- 
portant (see Fig. |2]l. Specifically, we can see from Fig. [T] 
that as the cm. energy is increased the ion-ion potential 
peak increases but the inner part of the barrier becomes 
narrower This is due to the fact that for high energies the 
system does not have enough time for rearrangements in 
the density to occur and the barrier approaches the frozen- 
density limit. However, at lower energies substantial den- 
sity rearrangements occur which modifies the inner part of 
the barrier. This modification is important for fusion cross- 
sections are deep sub-barrier energies. 



3.2 Deformed systems 

The collision of the ^"^Ni-i-'^'^Sn system represents a good 
example of a collision involving a deformed (oblate) nu- 
cleus, ^^Ni and a neutron rich nucleus. Fusion cross- 
sections for this system have been experimentally mea- 
sured |21| and initially a significant discrepancy was ob- 
served with standard coupled-channel calculations. We 
have used the DC-TDHF method to study this system ll22l 
|73|. The ion-ion potentials corresponding to two extreme 
orientations of the ''''Ni nucleus are shown in Fig.|3]as well 
as an empirical barrier used in barrier penetration calcula- 
tions in Ref. [21 1. Two important points are observed from 
this plot. The first is the strong dependence of the barrier 
height and location on the alignment of the deformed nu- 
cleus. We also see that the empirical barrier is very close to 
the equatorial orientation, which is closer to the assump- 
tion of spherical nuclei. The accuracy of our result with 
no parameters or normalization is impressive. The second 
point has to do with the meaning of sub-harrier, as seen 
from Fig. |3] while the experimental energies appear to be 
all sub-barrier with respect to the spherical barrier, two of 
them are above the barrier with respect to the fi - 90° 



potential barrier and the third faces a considerably nar- 
rower barrier. This explains the anomalous observation of 
enhanced fusion at these energies. Fig. |4] shows the ex- 
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Fig. 3. (Colour on-line) Potential barriers, ViR,P), obtained from 
DC-TDHF calculations for the ^''Ni+"^Sn system. Angle /? in- 
dicates different orientations of the deformed ^"'Ni nucleus in 
Ap = 10° intervals. Also shown are the experimental energies. 



perimental and theoretical fusion cross-sections calculated 
with different methods. The coupled-channel calculations 
are modified to include multiple neutron transfer As we 
see again the DC-TDHF results reproduce the fusion cross- 
sections reasonably well. 
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Fig. 4. (Colour on-line) Total fusion cross section as a func- 
tion of £c.m.- Shown are the experimental data (filled circles), the 
latest coupled-channel calculations |21 1 including neutron trans- 
fer (blue solid curve), and the DC-TDHF cross sections (dashed 
curve). 
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3.3 Superheavy systems 

Ion-ion interaction potentials calculated using DC-TDHF 
correspond to the configuration attained during a particu- 
lar TDHF collision. For light and medium mass systems as 
well as heavier systems for which fusion is the dominant 
reaction product, DC-TDHF gives the fusion barrier with 
an appreciable but relatively small energy dependence. On 
the other hand, for reactions leading to superheavy systems 
fusion is not the dominant channel at barrier top energies. 
Instead the system sticks in some dinuclear configuration 
with possible break-up after exchanging a few nucleons. 
For this reason the energy dependence of the DC-TDHF 
interaction barriers for these systems is not just due to the 
dynamical effects for the same final configuration but ac- 
tually represent different final configurations. For the same 
reasons calculations presented here can only address the 
capture cross-section for these systems since the long-time 
evolution to complete fusion or break-up is beyond the 
scope of TDHF due to the absence of quantum decay pro- 
cesses and transitions. 
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Fig. 5. (colour on-line) Potential barriers, V(R), for the 
'*'*Ca+^^^U system obtained from DC-TDHF calculations as a 
function of £c m. energy and for selected orientation angles yS of 
the ^^^'^U nucleus. Also, shown are the experimental cm. energies. 



As an example of superheavy formation from a hot- 
fusion reaction we have studied the "^^Ca-i-^^^U system 
||24l. Hartree-Fock (HF) calculations produce a spherical 
"^^Ca nucleus, whereas ^^^U has a large axial deformation. 
The large deformation of ^^^U is expected to strongly influ- 
ence the interaction barriers for this system. This is shown 
in Fig. |5] which shows the interaction barriers, V(R), cal- 
culated using the DC-TDHF method as a function of cm. 
energy and for three different orientations of the ^^^U nu- 
cleus. The alignment angle /3 is the angle between the 
symmetry axis of the ^^**U nucleus and the collision axis. 
Also shown in Fig. |5]is the point Coulomb potential cor- 
responding to this collision. The deviations from the point 
Coulomb potential at large R values are due to the defor- 
mation of the ^^^U nucleus. We first notice that the barriers 
corresponding to the polar orientation (J3 = 0°) of the ^^**U 



nucleus are much lower and peak at larger ion-ion separa- 
tion distance R. On the other hand, the barriers correspond- 
ing to the equatorial orientation of ^■'^U (J3 - 90°) are much 
higher and peak at smaller R values. For the intermediate 
values of /3 the barriers rise rapidly as we increase the ori- 
entation angle from /3 = 0°, as can be seen for jS = 45°. 
The rise in the barrier height as a function of increasing 
/3 values is not linear but seems to rise more rapidly for 
smaller/? values. We also see that for lower energies cen- 
tral collisions with polar orientation of ^^^U are the only 
orientations which result in the sticking of the two nuclei, 
while the equatorial orientations of ^^^U result in a deep- 
inelastic collision. Also, shown in Fig.|5]are the experimen- 
tal energies ||251 for this reaction. We observe that all of 
the experimental energies are above the barriers obtained 
for the polar alignment of the ^^^U nucleus. For the calcu- 
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Fig. 6. (colour on-line) Capture cross-sections for the ^^Ca+^^'^U 
system as a function of Ec.,,, energy (black circles). Also, shown 
are the experimental cross-sections (red squares) |25|. 



lation of the capture cross-section we need to average over 
all possible alignments of the ^-"^U nucleus as indicated by 
Eq. (fT4l l. Due to the relatively small charge of the "^^Ca nu- 
cleus the alignment probability dP/dQ of Eq. (fT4l i is in the 
range 0.46 - 0.52 and does not vary appreciably with en- 
ergy. In Fig.|6]we show the capture cross-sections for the 
"^^Ca-H^^^U system as a function of E^.m. energy (black cir- 
cles). Also, shown are the experimental cross-sections (red 
squares) |25|. One important fact to notice in the cross- 
section formula given in Eq. ( fT4l i is that the cross-section 
is multiplied by the sin(/3) factor, which renders the contri- 
bution originating from the lowest barriers at small values 
of yS to be very small. 

In Fig.|7]we also show the excitation energy E*(R) as 
a function of cm. energy and for two alignment angles 
(J3 = 0° and /? = 90° ) of the ^^^u nucleus. The excita- 
tion energy curves start at zero excitation when the two 
nuclei are far apart, which also provides a test for the nu- 
merical accuracy of the calculation. We note that the sys- 
tem is excited much earlier during the collision process for 
the polar aUgnment of the ^^**U nucleus and has a higher 
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Fig. 7. (colour on-line) Excitation energy, E'(R), in a central col- 
lision for various E^ ,n_ energies and for alignment angles angles 
/? = 0° and j8 = 90° of the ^''^U nucleus. 

excitation than the corresponding colUsion for the equato- 
rial orientation. Only two curves are shown for the equato- 
rial collision since at lower energies we have deep-inelastic 
collisions for this alignment. We note that the highest point 
reached for these excitation curves is chosen to be the point 
where the nuclei almost come to a stop inside the barrier, 
which corresponds to a nearly zero collective kinetic en- 
ergy. Since this is determined during the initial phase of 
the collision the dinuclear system is not in thermal equi- 
librium. However, the system essentially oscillates about 
this point. For energies for which the collision outcome is 
capture this would be the excitation energy at the capture 
point. 



4 Conclusions 

In summary, we have used the fully microscopic DC- 
TDHF method to obtain ion-ion interaction barriers for 
calculating fusion and capture cross-sections. The stan- 
dard, parameter free approach of DC-TDHF yields poten- 
tial barriers that can accurately reproduce the fusion cross- 
sections. The DC-TDHF approach has now been applied 
to study a number of systems with very promising results. 
This further elucidates the point that for the proper de- 
scription of fusion cross-sections dynamical effects such 
as neck formation and mass transfer must be included to 
modify the inner part of the potential barrier. Also interest- 
ing are the implications of the energy dependent barriers. 
These findings underscore the fact that additional modifi- 
cations needed for phenomenological methods to explain 
the fusion cross-sections may largely be due to the inade- 
quacy of the approximations made in treating the nuclear 
dynamics. 
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